Controlling spatiotemporal chaos in excitable media by local biphasic stimulation 
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Controlling spatiotemporal chaos in excitable media by applying low-amplitude perturbations locally is of immediate ap- 
plicability, e.g., in treating ventricular fibrillation, a fatal disturbance in the normal rhythmic functioning of the heart. We 
look at a mechanism of control by the local application of a series of biphasic pulses, i.e., involving both positive and negative 
stimulation. This results in faster recovery of the medium, making it possible to overdrive the chaos by generating waves with 
frequency higher than that possible with only positive pulses. This provides the simplest and most general understanding 
of the effectiveness of biphasic stimulation in controlling fibrillation and allows designing optimal waveshapes for controlling 
spatiotemporal chaos. 
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A characteristic feature of excitable media is the for- 
mation of spiral waves and their subsequent breakup into 
spatiotemporal chaos. Examples include catalysis of CO 
on Pt(llO) surface [1], cAMP waves during slime mold 
aggregation [2], etc. Another example of obvious im- 
portance is the propagation of waves of electrical exci- 
tation along the heart wall, initiating the muscular con- 
tractions that enable the heart to pump blood. In fact, 
spiral turbulence has been identified by several investi- 
gators as the underlying cause of certain arrhythmias, 
i.e., abnormal cardiac rhythms, including ventricular fib- 
rillation (VF) [3], a potentially fatal condition in which 
different regions of the heart are no longer activated co- 
herently. Current methods of defibrillatory treatment in- 
volve applying large electrical shocks to the entire heart 
in an attempt to drive it to the normal state. However, 
this is not only painful but also dangerous, as the re- 
sulting damage to heart tissue can form scars that act 
as substrates for future cardiac arrhythmias. Devising 
a low-amplitude control mechanism for spatiotemporal 
chaos in excitable media is therefore not only an excit- 
ing theoretical challenge but of potential significance for 
the treatment of VF. In this paper, we propose using 
low-amplitude biphasic pacing, i.e., applying a sequence 
of alternating positive and negative pulses locally, as a 
robust control method for such chaos. 

Most of the methods proposed for controlling spa- 
tiotemporal chaos in excitable media involve applying 
perturbations either globally or over a spatially extended 
system of control points covering a significant proportion 
of the entire system [4-7]. However, in most real situa- 
tions this may not be a feasible option. Further, in the 
specific context of controlling VF, a local control scheme 
has the advantage that it can be readily implemented 
with existing hardware of the Implantable Cardioverter- 
Defibrillator (ICD). This is a device implanted into pa- 
tients at high risk from VF that monitors the heart 
rhythm and applies electrical treatment, when necessary, 
through electrodes placed on the heart wall. A low- 



energy control method involving ICDs should therefore 
aim towards achieving control of spatiotemporal chaos by 
applying small perturbations from a few local sources. 

For most of the simulations reported in this paper we 
have used the modified Fitzhugh-Nagumo equations pro- 
posed by Panfilov as a model for ventricular activation 
[8] . For simplicity we assume an isotropic medium; in this 
case the model is defined by the two equations governing 
the excitability e and recovery g variables. 



de/dt = 
dg/dt = 



\/'e-f{e)-g, 
e{e,g){ke-g). 



(1) 



The function /(e), which specifies fast processes (e.g., 
the initiation of excitation) is piecewise linear: /(e) = 
Cie, for e < ei, /(e) = —€26 + a, for ei < e < e2, 
and /(e) = C3(e — 1), for e > e2. The function e(e,^), 
which determines the dynamics of the recovery variable, 
is e(e,^) = ei for e < e2, e(e,^) = 62 for e > e2, and 
e(e,^) = 63 for e < ei and g < gi. We use the physically 
appropriate parameter values given in Ref. [6]. 

We solve model (1) by using a forward-Euler inte- 
gration scheme. The system is discretized on a two- 
dimensional grid of points with spacing Sx = 0.5 dimen- 
sionless units. The standard five-point difference stencil 
is used for the 2-D Laplacian. The spatial grid consists 
of L X L points; in our studies we have used values of 
L upto 500. For the 1-D simulations we use a 3-point 
difference stencil for the Laplacian, with the spatial lat- 
tice consisting of L points. The time step is dt = 0.01 
dimensionless units. We did not observe any qualitative 
change in the results on decreasing the space and time 
steps by a factor of 2. On the edges of the simulation 
region we use no-fiux (Neumann) boundary conditions. 

To achieve control of spatiotemporal chaos, we locally 
apply a periodic perturbation, AF(27r/t), of amplitude 
A and frequency /. F can represent any periodic func- 
tion, e.g., a series of pulses having a fixed amplitude 
and duration, applied at periodic intervals defined by the 
stimulation frequency /. The control mechanism can be 
understood as a process of overdriving the chaos by a 
source of periodic excitation having a significantly higher 



frequency. As noted in Refs. [9,10], in a competition be- 
tween two sources of high frequency stimulation, the out- 
come is independent of the nature of the wave generation 
at either source, and is decided solely on the basis of their 
relative frequencies. This follows from the property of an 
excitable medium that waves annihilate when they collide 
with each other [11]. The lower frequency source is even- 
tually entrained by the other source and will no longer 
generate waves when the higher frequency source is with- 
drawn. Although we cannot speak of a single frequency 
source in the case of chaos, the relevant timescale is that 
of spiral waves which is limited by the refractory period 
of the medium, Tref, the time interval during which an 
excited cell cannot be stimulated until it has recovered its 
resting state properties. To achieve control, one must use 
a periodic source with frequency / > t~^\' This is almost 
impossible with purely excitatory stimuli as reported in 
Ref. [6] ; the effect of locally applying such perturbations 
is essentially limited by refractoriness to the immediate 
neighborhood of the stimulation point. 

A simple argument shows why a negative rectangu- 
lar pulse decreases the refractory period for the Panfilov 
model in the absence of the diffusion term. The stimu- 
lation vertically displaces the e-nullcline of Eq. (1) and 
therefore, the maximum value of g that can be attained 
is reduced. Consequently, the system will recover faster 
from the refractory state. To illustrate this, let us assume 
that the stimulation is applied when e > 62- Then, the 
dynamics reduces to e = —Cs{e — 1) — g^g = e2{ke — g). 
In this region of the (e,^)-plane, for sufficiently high ^, 
the trajectory will be along the e-nullcline, i.e., e :^ 0. If 
a pulse stimulation of amplitude A is initiated at t = 
(say), when e = e(0),^ = ^(0), at a subsequent time t. 
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where, a = e2k{l + ■^), 6 = e2[l + (k/Cs)]- The negative 
stimulation has to be kept on till the system crosses into 
the region where e < 0, after which no further increase 
of g can occur, as dictated by the dynamics of Eq. (1). 
Now, the time required by the system to enter this re- 
gion is ^ln^^y^3^, where (j) = Cs{l — ^2) +^- Therefore, 
this time is reduced when A < and contributes to the 
decrease of the refractory period. Note that, the above 
discussion also indicates that a rectangular pulse will be 
more effective than a gradually increasing waveform, e.g., 
a sinusoidal wave (as used in [12]), provided the energy 
of stimulation is same in both cases, as the former allows 
a much smaller maximum value of g. Therefore, phase 
plane analysis of the response to negative stimulation al- 
lows us to design waveshapes for maximum efficiency in 
controlling spiral turbulence. 

To understand how negative stimulation affects the re- 
sponse behavior of the spatially extended system, we first 
look at a one-dimensional system. Fig. 1 shows the re- 
lation between the stimulation frequency / and effective 
frequency /e//, measured by applying a series of pulses 



at one site and then recording the number of pulses that 
reach another site located at a distance without being 
blocked by any refractory region. Depending on the rel- 
ative value of / and r^e/, we observe instances of n : m 
response, i.e., m responses evoked by n stimuli. From 
the resulting effective frequencies /e//, we can see that 
for purely excitatory stimulation, the relative refractory 
period can be reduced by increasing the amplitude A. 
However, this reduction is far more pronounced when a 
negative stimulation is applied between every pair of pos- 
itive pulses. The inset in Fig. 1 shows that there is an 
optimal time interval between applying the positive and 
negative pulses that decreases the refractory period by 
as much as 50%. The highest effective frequencies corre- 
spond to a stimulation frequency in the range 0.1 — 0.25, 
agreeing with the optimal time period of 2-5 time units 
between positive and negative stimulation. 

A response diagram similar to the one-dimensional 
case is also seen for stimulation in a two-dimensional 
medium (Fig. 2). A small region consisting of n x n 
points at the center of the simulation domain is chosen 
as the stimulation point. For the simulations reported 
here n = 6; for a smaller n, one requires a perturbation 
of larger amplitude to achieve a similar response. To 
understand control in two dimensions, we find out the 
characteristic time scale of spatiotemporal chaotic activ- 
ity by obtaining its power spectral density (Fig. 2, inset). 
We observe a peak at a frequency fc — 0.0427. As seen 
in Fig. 2, there are ranges of stimulation frequencies that 
give rise to effective frequencies higher than this value. 
As a result, the periodic waves emerging from the stim- 
ulation point will gradually impose control over the re- 
gions exhibiting chaos. If / is only slightly higher than fc 
control takes very long; if it is too high the waves suffer 
conduction block at inhomogeneities produced by chaotic 
activity that reduces the effective frequency, and control 
fails. Note that, at lower frequencies the range of stimu- 
lation frequencies for which /g// > /c, is smaller than at 
higher frequencies. We also compare the performance of 
sinusoidal waves with rectangular pulses, adjusting the 
amplitudes so that they have the same energy. The for- 
mer is much less effective than the latter at lower stimu- 
lation frequencies, which is the preferred operating region 
for the control method. 

The effectiveness of overdrive control is limited by the 
size of the system sought to be controlled. As shown in 
Fig. 3, away from the control site, the generated waves 
are blocked by refractory regions, with the probability of 
block increasing as a function of distance from the site of 
stimulation [13]. To see whether the control method is 
effective in reasonably large systems, we used it to termi- 
nate chaos in the two-dimensional Panfilov model, with 
L = 500 [14]. Fig. 4 shows a sequence of images illus- 
trating the evolution of chaos control when a sequence 
of biphasic rectangular pulses are applied at the center. 
The time necessary to achieve the controlled state, when 



the waves from the stimulation point pervade the entire 
system, depends shghtly on the initial state of the sys- 
tem when the control is switched on. Not surprisingly, 
we find that the stimulation frequency used to impose 
control in Fig. 4 belongs to a range for which /g// > fc 

Although most of the simulations were performed with 
the Panfilov model , the arguments involving phase plane 
analysis apply in general to excitable media having a 
cubic-type nonlinearity. To ensure that our explanation 
is not sensitively model dependent we obtained similar 
stimulation response diagrams for the Karma model [15]. 

Some local control schemes envisage stimulating at 
special locations, e.g., close to the tip of the spiral wave, 
thereby driving the spiral wave towards the edges of the 
system where they are absorbed [16]. However, aside 
from the fact that spatiotemporal chaos involves a large 
number of coexisting spirals, in a practical situation it 
may not be possible to have a choice regarding the loca- 
tion of the stimulation point. We should therefore look 
for a robust control method which is not critically sensi- 
tive to the position of the control point in the medium. 
There have been some proposals to use periodic stimula- 
tion for controlling spatiotemporal chaos. For example, 
recently Zhang et al [12] have controlled some excitable 
media models by applying sinusoidal stimulation at the 
center of the simulation domain. Looking in detail into 
the mechanism of this type of control, we have come to 
the conclusion that the key feature is the alternation be- 
tween positive and negative stimulation , i.e., biphasic 
pacing, and it is, therefore, a special case of the general 
scheme presented here. 

Previous explanations of why biphasic stimulation is 
better than purely excitatory stimulation (that use only 
positive pulses), have concentrated on the response to 
very large amplitude electrical shocks typically used in 
conventional defibrillation [17,18] and have involved de- 
tails of cardiac cell ion channels [19]. To the best of 
our knowledge the present paper gives the simplest and 
most general picture for understanding the efficacy of the 
biphasic scheme using very low amplitude perturbation, 
as it does not depend on the details of ion channels re- 
sponsible for cellular excitation. 

There are some limitations to achieving control over a 
large spatial domain in an excitable medium by pacing at 
a particular point. Under some parameter regimes, the 
circular waves propagating from this point may them- 
selves become unstable and undergo conduction block at 
a distance from the origin (similar to the process out- 
lined in Ref. [20]). In addition, the control requires a 
slightly higher amplitude and has to be kept on for peri- 
ods much longer than spatially extended control methods 
[6]. However, these drawbacks may be overcome if we use 
multiple stimulation points arranged so that their regions 
of inffuence cover the entire simulation domain. 

In conclusion, we have proposed a simple explanation 
of the efficacy of low-amplitude biphasic stimulation in 



controlling spatiotemporal chaos in excitable media (e.g., 
VF). It is based on the competition between the fre- 
quency of the applied stimulation, and the effective fre- 
quency (obtained from the dominant timescale) of chaos. 
The former can be increased relative to the latter only 
by decreasing the refractory period which is achieved by 
a negative stimulus prior to applying the excitatory posi- 
tive stimulation. Our analysis makes it possible to design 
pacing waveforms for maximum efficiency in controlling 
chaos. 
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FIG. 1. Stimulation response diagram for one-dimensional 
Panfilov model (L = 40) for different stimulation frequencies 
/. The dotted and broken curves represent purely excitatory 
pulses of amplitude A = 5, pulse duration r = 0.05/~^, and 
A = 10, r = 0.1/~^, respectively, while the solid curve repre- 
sents biphasic pulses of amplitude A = 10 and pulse duration 
T = 0.1/~^ (as shown in the bottom left corner). Note that, 
the highest effective frequencies feff for the the three cases 
are very different. The ratio / : feff is shown for the first 
four peaks. The inset shows the decrease in refractory period 
in absence of the diffusion term when a negative stimulation 
is applied at different times (tneg) after the initial excitation. 
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FIG. 2. Stimulation response diagram for two-dimensional 
Panfilov model (L = 26) showing relative performance of dif- 
ferent waveforms. The dash- dotted line represents a sinus 
wave {A = 6) and the soUd curve represents a wave of bipha- 
sic rectangular pulses {A = 18.9), as shown in the bottom left 
corner, such that they have the same total energy. The inset 
shows the power spectra of spatiotemporal chaos in the 2-D 
Panfilov model (L = 500). The e variable was recorded for 
3300 time units and the resulting power spectral density was 
averaged over 32 points. The peak occurs at the characteristic 
frequency fc — 0.0427 which is indicated in the main figure 
by the broken line. 
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FIG. 3. Distance dependence of stimulus response for dif- 
ferent stimulation frequencies / in the two-dimensional Pan- 
filov model. Biphasic rectangular pulses {A = 18.9) having 
duration r = 0.1/~^ are applied, which elicit a response hav- 
ing an effective frequency feff at a particular location. The 
first three cells (x = 1,2,3) are within the region subject to 
direct stimulation. The shaded regions represent different re- 
sponse ratios / : feff- Integral multiples of the characteristic 
frequency fc are indicated on the /-axis. 
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FIG. 4. Control of spatiotemporal chaos in the two- di- 
mensional Panfilov model (L = 500) by applying biphasic 
pulses with amplitude A = 18.9 and frequency / = 0.13 at 
the center of the simulation domain. The pulse shape is rect- 
angular, having a duration of ^ 0.77 time units. Snapshots 
are shown for (top left) t = 0, (top right) t = 1000, (bottom 
left) t = 2700 and (bottom right) t = 3800 time units. The 
excitation wavefronts are shown in white, black marks the re- 
covered regions ready to be excited, while the shaded regions 
indicate different stages of refractoriness. 



